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In this paper, an optimal artificial neural network (ANN) controller for load 
frequency control (LFC) of a four-area interconnected power system with 
non-linearity is presented. A feed forward neural network with multi-layers 
and Bayesian regularization backpropagation (BRB) training function is 
used. This controller is designed on the basis of optimal control theory to 
overcome the problem of load frequency control as load changes in the 
power system. The system comprised of transfer function models of two 
thermal units, one nuclear unit and one hydro unit. The controller model is 
developed by considering generation rate constraint (GRC) of different units 
as a non-linearity. The typical system parameters obtained from IEEE press 
power engineering series and EPRI books. The robustness, effectiveness, 
and performance of the proposed optimal ANN controller for a step load 
change and random load change in the system is simulated through using 
MATLAB-Simulink. The time response characteristics are compared with 


that obtained from the proportional, integral and derivative (PID) controller 
and non-linear autoregressive-moving average (NARMA-L2) controller. The 
results show that the algorithm developed for proposed controller has a 
superiority in accuracy as compared to other two controllers. 
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1. INTRODUCTION 

Automatic load frequency control is the main area of concern in the operation of an interconnected 
power system. As load on the system changes, the frequency changes. In order to balance the megawatt 
(MW) generation with load demand, it is necessary to control the synchronous generators power output and 
frequency. This control is automatic, as it maintains the frequency at base value and power flows via tie-lines 
within scheduled values for a perturbation in the load [1]—[3]. The load frequency control monitors the area 
control error (ACE), which is the net real power interchange between the control areas plus the frequency 
deviation multiplied with a frequency bias. To reduce ACE, closer to zero, change the position of speed 
changer of the generator governors within the control area by means of load frequency control (LFC). 

In the design of load frequency controllers, conventional control techniques have been developed, 
which gives slower responses. The developments in advanced technology, artificial intelligent (AT) based 
techniques such as neural networks, fuzzy logic, genetic algorithm (GA) and particle swarm optimization 
(PSO) have been using and these techniques overcome the disadvantages of conventional controllers and 
increase the speed of response. When the load changes suddenly, the primary automatic load frequency 
control can be obtained by the action of speed governor in the prime movers. A supplementary or secondary 
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automatic load frequency control action is used in the proposed system to bring the change in frequency and 
change in tie-line power to zero with the help of advanced controllers. The past studies in the literature 
explain the LFC problem with proportional-integral-derivative (PID) and fuzzy logic controllers, but less 
with optimal artificial neural network (ANN) controller and nonlinearity. 

Prasad and Ansari [4], employed a three-layer ANN observer-based control strategy used in a two 
similar area interconnected power system with generation rate constraint (GRC) and governor dead band 
(GDB). This system is simulated for random unmatched disturbance estimation and its rejection. Qian and 
Fan [5], implemented a three-layer radial basis function (RBF) neural network for load frequency control of a 
two-area power system with GRC and wind turbine model. The control scheme is designed on the basis of 
terminal sliding mode control. Bhatia et al. [6] proposed a three-layer neural network-based NARMA-L2 
controller for a three similar area power system with GRC. Only frequency deviation is discussed. Chettibi 
et al. [7] proposed and implemented a technique for forecast of grid voltage frequency in short time based on 
ANN models and deep recurrent neural networks. This can be used in an advanced control scheme and 
monitoring distributed generators for frequency and voltage variations. The performance of these networks 
was assessed in terms of root mean square error (RMSE) that was lies between 0.002 Hz and 0.01 Hz for 
sapling interval of 0.1 s and 1 s respectively. 

Alzaareer et al. [8] proposed an ANN based NARMA-L2 controller model for a three-area 
interconnected system without GRC. This controller is compared with PI and PID controller for load 
frequency control in the system. Prakash and Sinha [9], proposed a hybrid neuro fuzzy (HNF) controller in a 
four-area power system without GRC. This controller performance is compared with fuzzy logic, three-layer 
ANN and PID controllers for 1% change in load. Peak overshoot and settling time values of -0.055 pu 
(-2.75 Hz) and 40 s respectively, obtained with ANN controller. Kumari et al. [10], proposed an ANN-PID 
control technique for a two area non-reheat thermal plant power system without GRC. The controller 
performance is tested with 1096 step load perturbation and different error values are measured. Prakash and 
Sinha [11], proposed an ANN and adaptive neuro fuzzy inference system (ANFIS) for a six-area power 
system composed of hydro, thermal, gas, diesel, and nuclear plants without GRC. The controller performance 
is tested with 1% step load perturbation. Mucka et al. [12], employed a three-layer neural network-based 
NARMA-L2 controller for a four-area power system without GRC. The system is simulated with 2% change 
in load at frequency 50 Hz and its response has more settling time and undershoot. 

The above research work [4]-[7], employs ANN based controller for two-area and three-area 
interconnected power system with GRC. Only two input variables and first order governor-turbine transfer 
function models are considered and the work [8]-[12], even though propose ANN based controllers with 
Levenberg-Marquardt learning function but does not provide information on the number of neurons in the 
hidden layer(s) and considered only first order non-reheat turbines without GRC. Hence, the present work 
proposes an optimal ANN controller with BRB training function and is designed based on state space model 
for load frequency control in a four-area power system comprises reheat tandem compound turbines with 
GRC and IEEE standard parameters are chosen within the operating constraints of system components. 


2. MODELLING OF THE SYSTEM 

The main components of each area are governor, turbine, generator, and load. The dynamic models 
of governors, tandem compound steam turbines and hydro turbines were presented in [13], [14]. For thermal, 
nuclear, and hydro power plants, the transfer function models of a governor or hydraulic valve actuator are 
obtained from the basic Watt's governor operation. The thermal plant governor and turbine block diagram 
with fraction of power generated by high pressure (HP), intermediate pressure (IP), and low pressure (LP) 
sections is shown in Figure 1. Figure 2 shows the block diagram of nuclear plant governor and turbine with 
fraction of power generated by very high pressure (VHP), high pressure (HP), and low pressure (LP) 
sections. Figure 3 shows the block diagram of hydro plant governor and turbine. 

The generation rate constraint is the limitation on the rate of change in the real power generation due 
to physical limitations of turbine. The existence of GRC [5], [7], [15] has an adverse effect on system 
stability. It should be considered for LFC problem as a non-linear model shown in Figure 4. The GRC values 
are taken into account by adding limiters to the turbines. The GRC values for thermal and nuclear plants are 
+0.005 pu. MW.s~* and that for hydro plant is +0.045 pu. MW.s7~* and —0.06 pu. MW.s^!. The transfer 
function model of synchronous generator and load is obtained by rotor dynamics, swing equation and overall 
frequency dependent characteristic of a composite load. Synchronous generator-load transfer function model 
[1], [2] in standard first order form is obtained as (1). 


Gs, (s) = ———— = —PB4. fori = 1,2,3,4 (1) 


2SHj+Dj  1*STpsi 
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Figure 3. Hydro plant governor and turbine model 


Figure 4. Generation rate constraint model 


The output of generator-load model is the change in frequency or the frequency deviation Af;(s) 
due to change in load AP,;(s). In normal operation, the change in tie-line power is obtained from 
synchronizing torque coefficient (T) using (2). ACE [9] is the input signal to controller for each power 
system area and is calculated using (3). 


2nT 


ZT (Afi (s) — Af; (s)) for i 212,34 (2) 


S 


AP (s) = 


The objective function (OF) determines the system dynamics and satisfy criterion such as fast response with 
minimized undershoot and steady state error. Thus, integral of time weighted absolute error (ITAE) is used as 
OF [15] and is calculated as (4). 


ITAE; = f” t|AP; + B;Af;| dt (4) 


0 


By connecting the block diagrams of governor, turbine, generation rate constraint, generator and load models 
of respective areas and interconnecting these areas via tie-line model gives the complete block diagram of a 
four-area interconnected power system as shown in Figure 5. 
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Figure 5. Four area interconnected power system 


3. PID AND NARMA-L2 CONTROLLERS 
3.1. PID controller 

The PID controllers are conventional controllers used when the system requires improvement under 
steady-state and transient conditions. These controllers design is simple and inexpensive. The Ziegler- 
Nichols method proposed in [16], [17], is employed to determine the tuned gain values of proportional (Kp), 
integral (K;) and derivative (K4). 


32. NARMA-L2 controller 

The non-linear autoregressive moving average controller is the most effective in the non-linear 
control systems. It is referred to as NARMA-L2 control when the plant model can be approximated by a 
particular form [18]-[21]. The dynamic responses of the area frequency and tie-line power flows are obtained 
using this controller in the power systems [22]. Its main function is to transform non-linear system dynamics 
into linear dynamics by cancelling the non-linearities. 


4. RESEARCH METHOD-OPTIMAL ANN CONTROLLER 

In the design of load frequency optimal controller, an artificial neural network (ANN) is to be 
trained. The flow chart of neural network training process is shown in Figure 6. The training process is 
divided into three main sections, which are pre-training steps, training the network, and post-training 
analysis. 
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Figure 6. Flow chart of neural network training process 


4.1. Data collection and preprocessing 
The optimal controller is designed for the power system using state space model [23] with 32 state 
variables and 4 control output variables. The aim of this controller is to obtain a control law u(x,t) for 
minimizing the performance index. Formulation of the state space model is achieved by writing differential 
equations [13], [24] describing each individual block of four area power system in terms of state variables. 
These variables are output of various blocks represent the change in mechanical power, electrical 
power and frequency and are defined as: 


- State variables: 
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X = Afi, X2 = AP, X3 = APcgi; X4 = APggi X5 = Acur Xo = APy,, X7 = Af, Xg = APmo2, 


X9 = APco2, X10 = APg2, X11 = APp,, X12 = APcu2, X13 = APy2, X14 = Afo, X15 = AP, 


X16 = APco3, X17 = APgga, Xig = APcys, X19 = APy3, X20 = Afa X21 = APma, X22 = APur, 


X23 = APyp, X24 = APgg, X25 = AP 2 + AP + AP y4, X26 = AP5; + APz3 + APo4, 


X5; = AP3, + AP3; + AP34, X28 = AP}, + AP + AP X09 = f ace, dt, X39 = | ace. dt 


X31 = | ace, dt, x35 = face, dt 


Control inputs: 14, u5, u3 and u4, disturbance inputs: di = APj4, d; = AP,2,d3 = AP;3 and d; = AP;4 


- State equations: 


For block 1: x 


For block 2: x 


For block 3: 


For block 4: x 


For block 5: 


For block 6: 


1 Kps1 Kps1 s1 
X4 == X4 + P X2 LIE X25 p 
Tps1 Tps1 Tps1 T, 
X2 = X3 
J 1 Kt Kti Kti Kth Kei Kth Kth 
ty =- x + ( +=- =) x, + + —# — Xs d — xe 
Ttc Tte Tte Ter Ttc Ter Ttt Ttt 
1 d 1 
X4 = CU. Xs an: Xs = — — Xs 7, *6 
tr tt 
1 1 1 
Xg =— = Xg + Uy 
Ri Ttg Ttg Ttg 
4 1 Kps2 Kps2 Kps2 
X7 Im X7 + 4 Xg E X26 -2 d 
Tps2 Tps2 Tps2 T, 
Xg E Xo 
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(7) 
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(9) 


(10) 
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: 1 K K K K K K K 
For block 7: X9 = — — xg + zx, + ( nh ut) X11 + ( ny ne 4 ER TT xis 
Tnc Tnc Tnc Tnrı Tne Tnt Tnri Tnt 
T 1 1 AME 1 1 
Xio = ———JX49 t L—X4.X41— ———3X41 + —— X42 and 
Tnr2 Tnr2 Tnr1 Tnri 
X12 — RT. +7 X13 (11) 
: 1 1 
For block 8: x44 — — X; X43 t —u5 (12) 
R2 Tng Tng Tng 
; 1 K K K 
For block 9: X44, = — ——2x44 + 22 xis — P2 x — 22d (13) 
Tpsa Tps3 Tps3 Tps3 
For block 10: xis = x16 (14) 
. 1 K Kti Kg K Ky K K 
For block 11: 3,5 = — x, + ($E + E — St) x + (Eth + Rt LR) ag + xe 
Ttc Ttc Tte Ter Ttc Ttr Ttt Ttt 
. 1 1 : 1 
X1; = —— X17 t —Xıg and žig = — — X48 + — X19 (15) 
Ttr Ttr Ttt Ttt 
: 1 1 1 
For block 12: X49 = — X14 — = X19 t — us (16) 
R3 Ttg Ttg Ttg 
; 1 K K K 
For block 13: X29 = — ——x?9 + 2x44 — P xs — Pd, (17) 
Tps4 Tps4 Tps4 Tps4 
For block 14: x3, = x22 (18) 
; 1 1 
For block 15: x5; = — | ————— 4s J X29 — 4——— X22 + 
0.5 R4 Thg(1+ ait) 0.5 Thw 


1- Inr 
1 1 Thg 1 
"TEM + os r, [14 Ft) X23 — WE 14 REY X24 — osr, [afre u, (19) 
w o Thr(1+ Tee) 0. Ty (14 c) 0. Thg(1+ et) 


1 1 
For block 16: 335 > cm X20 — (ao) X23 + 
44hg Rh. hr Rn. 


T; 
1- hr 
T, 1 
T X24 + Rnt d 
qu (1 + zi) Tng (1 + r9 
em ed e: 1 
X24 — RaThg X20 Tng X24 + Tng ug (20) 


For tie-lines: X95 = 2m T (3x1 — x; — X14 — X20), X26 = 2t T (3x7 — x4 — X14 — X20) 
X5; = 2n T(3x44 — X4 — X; — X20), *28 = 2T T (3X20 — x4 — x; — x14) Q1) 


For controller inputs: X29 = By x, + X25, X30 = Bx; + X26> X31 = B4x44 + X27, 
X32 = B4X29 + X28 (22) 


Then, the state equation in matrix form: 

ž = Ax + Bu + Fd (23) 
Output equation: 

y =Cx (24) 
where the matrix A (32x32) is a coefficient matrix of all the state variables, the matrix B (32x4) is a 
coefficient matrix of all the control variables, the matrix F (32x4) is a coefficient matrix of all the disturbance 
variables, the matrix C (1x32) is a coefficient matrix of output variables, x = [x4, X2, ...,X32]"= state vector, 


u = [u ..u4]T 7 control vector and d = [d, ...d,]"=disturbance vector. 
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The optimal control inputs vector, u = —Kx is obtained by a linear combination of all states, where 
K is the feedback gain matrix. MATLAB code is used to obtain the matrix K by solving of the reduced matrix 
Riccati equation [6], [9], [20] given by (25): 


ATS + SA — SB[R-1B" 5] +Q =0 Q5) 


where R^! BT'S = K and matrix S is a real, positive definite and symmetric. The matrices Q and R are 
determined on the basis of three considerations: the excursions of ACE's, f ACE s dt and control inputs 
Uz ...u4 about steady values are minimized. These can be recognized as symmetric matrices to minimize 
performance index in quadratic form, given by (26) and (27). 


PI =~ 


=> fo GT Qx * xTRu)dt (26) 
— (B1x1)? + 2B1x1X25 + (x25)? + (B5x7)? + 2B2X7xX26 + (x26)? + (B3X14)? 

PI = zh +2B3x14%27 + (x27)? + (B4x29)* + 2B4X29X28 + (X28)? + (x29)? + dt (27) 
(xao)? + (x31)? + (x32)? + (u)? + (uz)? + (ua)? + (u4)? 


The discretized system state equations and optimal control inputs vector are used to collect/generate 
the training data for different values of step load change. Since the time of study and sampling have been 
chosen as 90 s and 0.005 s respectively, a total of 9000 samples are collected for each variable for a step load 
change simultaneously in all the four areas. All such variables form one data set, comprises of 40 variables 
(X4, X2, ..., X32, dy ... d4, Uy .. U4 ). Two data sets for each load disturbances have been collected. 


4.2. Selecting the neural network architecture 

A multilayer feedforward neural network architecture [25] shown in Figure 7 is employed for LFC 
in a non-linear four-area interconnected power system. The input scalar vector p is represented by a vertical 
bar with R inputs. There are 36 input nodes (R=36) corresponding to two input parameters which are 32 input 
nodes, each corresponding to 32 state variables x4, x2, ..., X32, and another 4 input nodes for load disturbances 
or perturbations (d, ...d4) in the system. The two hidden layers with hyperbolic tangent sigmoid transfer 
function is used with S neurons to verify the dependency of state variables with the load perturbations and the 
repeatability of convergence. Hidden layer 1 has St = 20 neurons, hidden layer 2 has S? = 10 neurons, and 
the output layer has S? — 4 neurons with linear transfer function are used in the network. The outputs of 
hidden layers 1 and 2 are the inputs for hidden layer 2 and output layer respectively. The vectors 
nl,n? and a? represent the net inputs, and at, a? and a?represent the outputs of hidden layers 1, 2 and 
output layer, respectively. For hyperbolic tangent sigmoid transfer function, input/output relations are given 
by (28). 


2..€ =e 
a aie cae and a*= aca (28) 
Inputs Hidden Layer 1 (Tan-Sigmoid) Hidden Layer 2 (Tan-Sigmoid) Outputs (Linear) 
Figure 7. Neural network architecture 
For linear transfer function, input/output relation is given by, a? — n? (29) 
The outputs of hidden layer 1: za! = tan sig(Wpxb?) (30) 
The outputs of hidden layer 2: a? = tan sig(W?a!xb?) (31) 
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The outputs of output layer: a? = purelin(W? a*xb?) (32) 


where, W1,W? and W? represent the weight matrices and bt, b? and b? represent bias vectors of hidden 
layers 1, 2 and output layer respectively. The output of network a?is a (4X1) vector represents the 
target/control signals (u, ... u4) given to power system corresponding to disturbances (d, ... d4). 


4.3. Neural network training 

The network training is based on the data collected from the optimal controller for different 
perturbations or step load changes. As the inputs are applied to the network, the outputs are compared to the 
target values, and the supervised learning rule is used. Before network training, initialize the weights and 
biases using the method of Widrow and Nguyen. In this method, set row i of Wt, ;W?!, to have a random 
direction and a magnitude of | WW! | =0.7(S 1/ R and set b; to a uniform random value between -| w!| 


and | wt]. Then, Bayesian regularization backpropagation is used to train the network for 100 epochs 
because it is very effective algorithm for training multilayer networks and generalized without the need for 
the validation data set. According to Levenberg-Marquardt optimization, weights and bias values are 
updated, as it minimizes a combination of squared errors and weights. The performance or mean squared 
error (MSE) goal is set to a very low value of 1 x 107t°. Marquardt adjustment parameter and minimum 
performance gradient are 0.005 and 1 x 10~’respectively. The training is carried out till the completion of 
100 epochs or MSE reaches the desired limit. The network was trained with different number of neurons in 
the hidden layers. The best MSE value (9.417 x 1071) is obtained at epoch 45, as shown in the Figure 8. 
This algorithm computes the effective number of parameters (y = 220) that are being used by the network. 
The values of MSE and y indicate that the network with St = 20 and S? = 10 is satisfactory. 


Best Training Performance is 9.417e-11 at epoch 45 


Mean Squared Error (mse) 


0 5 10 15 20 25 30 35 40 45 
45 Epochs 


Figure 8. Neural network training error curve 


5. RESULTS AND DISCUSSION 

The parameter values of system components are given in Table 1 at base frequency of 50 Hz. The 
PID controller is tuned by Z-N method and its parameters are shown in Table 2. MATLAB-Simulink is used 
to perform the simulation of a four-area interconnected power system with three types of controllers. Each of 
the parameters obtained from the range given in IEEE press power engineering series and EPRI books [26], 
[27] based on the design of the electrical components for best performance. 


Table 1. Power system parameter values 


Parameters P. P, Prie fo Tie ó R 
All Areas 2000 MW 1000 MW 200 MW 50 Hz 0.0866 30° 2.5 Hz/pu MW 
Thermal Tig Tet Ter Tie Kin Ky Ki D B H Kps Tps 
Plant 0.2s 0.3 s 7s 0.4s 0.3 0.4 0.3 0.01 0.41 5 100 20s 
Nuclear Tag Tat Tarı Tarz Tac Kuv Ku Kan = = a E 
Plant 0.25 0.3s 7s 78 0.4s 0.22 0.56 0.22 0.01 0.41 5 100 20s 
Hydro Thg Thw Thr Rtn Ra 5 = = s = Š i 
Plant 10s Is 5s 0.2875 0.05 - - 0.015 0.415 4 66.6667 10.6667s 
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Table 2. PID controller parameters 


Plant Ke Po K, Ki Ka 
Thermal 0.286 12.289 0.1716 0.0279 0.2636 
Nuclear 0.1814 19.137 0.1088 0.0114 0.2604 

Hydro 0.1119 16.885 0.0671 0.0080 0.1417 


MATLAB code is written to obtain optimal controller gain values, to train and test the proposed 
optimal ANN controller for load frequency control. A step change in load power AP, (steps of 1% up to 5%) 
in each area is applied. The time domain characteristics-settling time (t,) and undershoot (Mp), and errors of 
all four areas are measured and tabulated in Table 3. These specifications are measured using MATLAB 
function. Figure 9 shows the comparison of frequency deviations in area-1 and area-4 with different 
controllers under +2% step load change in each area. From Figures 9(a) and 9(b), the change in settling time 
and maximum undershoot values are measured and tabulated in Table 3. As the load on the system increases 
suddenly, its frequency decreases at that moment. The maximum undershoot decreases, oscillates, and settle 
to zero steady state value quickly due to optimal ANN controller action compared to other two controllers. 
This indicates that the system is stable even with +5% change in load. 

It is evident from the Table 3 that the proposed optimal ANN controller gives stable responses with 
a very minimum steady state error, lesser undershoot, lower settling time and the very smaller values of 
MSE. Figure 10 shows the frequency deviation with optimal ANN controller under equal and unequal load 
increase. Figure 10(a) shows for equal load (APj, = APj; = AP,3 = APj4 = 4%) in each area. Also, from 
the Table 3, settling time is 21.8343 s and undershoot is -0.2572 Hz with a MSE value of 1 x 1071?. These 
values are smaller compared to that with PID and NARMA-L2 controllers [19], [20] for the same change in 
load. Figure 10(b) shows the frequency deviation with optimal ANN controller under unequal load increase 
(AP, = 1%, AP;; = 2%, AP, = 396, AP, = 4%) in each area. Under this condition, proposed optimal 
ANN controller gives good dynamic response with zero steady state error. 

The load and frequency deviations in Figure 11 shows the robustness of the proposed optimal ANN 
controller. Figure 11(a) shows a random load pattern and is more realistic in real power systems. All four 
areas are encounter this type of load variations [15]. Under this condition, frequency deviation in all areas is 
shown in Figure 11(b). It reveals that, optimal ANN controller successfully tracks the load pattern and balance 
generation with load effectively with constant frequency. Only frequency deviation occurs for the equal change 
in percentage of load in each area, whereas the algebraic sum of change in tie-line power flow is zero. 

With optimal ANN controller in the four-area interconnected system, -0.0713 Hz and 23 s are the 
maximum values of undershoot and settling time for +1% change in load, respectively. For +5% load change, 
the peak undershoot is -0.3566 Hz and settling time is 22.9970 s. It is seen from the responses with 296 
increase in load causes a minimum undershoot of -0.1286 Hz and minimum settling time of 21.864 s. It is 
observed that, the settling time is constant as the step load increases from 1% to 5%. The magnitude of 
frequency deviation increases with load, but this increase is very small. On the other hand, as load decreases, 
the frequency deviation increases with the same settling time. The ITAE values measured with PID and 
NARMA-L2 controllers [20] are high compared to mean squared error values measured with optimal ANN 
controller. The time response specification values obtained are smaller compared to their values in literatures 
[28]-[32]. 


Table 3. Comparative study of settling time, maximum undershoot and error 


AP, Controllers Settling time (s) and error Maximum undershoot (Hz) and error in frequency 
in frequency deviation deviation 
Area-|  Area-2  Area-3 Area-4 ITAE/MSE  Area-l Area-2  Area-3 Area-4  ITAE/MSE 
PID 37.5065 37.1265 37.5065 36.8291 0.1128 -0.0741 -0.0752 -0.0741 -0.0749 0.1128 
+1% NARMA-L2 15.1145 15.0993 15.1078 26.6548 le-3 -0.0670 -0.0681 -0.0670 -0.0705 le-3 
ANN 21.8390 23.0015 21.8396 21.9799 le-10 -0.0643 -0.0656 -0.0643 -0.0713 le-10 
PID 31.5815 31.5713 31.5815 31.4318 0.3229 -0.2334 -0.2341 -0.2334 -0.2358 0.3229 
+2% NARMA-L2 18.1999 18.4093 18.2013 18.9832 le-3 -0.2297 -0.2297 -0.2297 -0.2336 le-3 
ANN 21.8640 23.0100 21.8640 21.9904 le-10 -0.1286 -0.1312 -0.1286 -0.1427 le-10 
PID 25.8792 25.8051 25.8792 25.8254 0.8842 -0.4649 -0.4649 -0.4649 -0.4714 0.8842 
+3% NARMA-L2 23.0941 23.1631 23.0240 22.8366 le-3 -0.4582 -0.4582 -0.4582 -0.4657 le-3 
ANN 21.8428 22.9973 21.8439 21.9854 le-10 -0.1929 -0.1968 -0.1929 -0.2140 le-10 
PID 48.4521 49.0057 48.4521 47.7849 1.9070 -0.7439 -0.7438 -0.7439 -0.7502 1.9070 
+4% NARMA-L2 31.4917 31.7661 31.4906 32.1603 le-3 -0.7284 -0.7284 -0.7284 -0.7319 le-3 
ANN 21.8347 22.9910 21.8343 21.9830 le-10 -0.2572 -0.2624 -0.2572 -0.2853 le-10 
PID 53.3577 53.7863 53.3577 52.8282 3.5370 -1.0526 -1.0526 -1.0526 -1.0549 3.5370 
+5% NARMA-L2 33.3008 33.5357 33.4160 33.9808 le-3 -1.0146 -1.0146 -1.0146 -1.0222 le-3 
ANN 21.8419 22.9970 21.8414 21.9850 le-10 -0.3215 -0.3280 -0.3215  -0.3566 le-10 
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Figure 9. Comparison of frequency deviations in (a) area-1 and (b) area-4 with 296 load increase 
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Figure 10. Frequency deviation with optimal ANN controller under (a) +4% load change in each area and 
(b) unequal load increase in the areas 


Df- Area-1 


== = Df-Area-2 


EE Df- Area-3 


-— Df- Area-4 


Random load change ( pu ) 
S 
s 
22] 

Frequency deviation ( Hz ) 


-0.015 


-0.025 
0 


20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180 
Time (s) Time (s) 
(a) (b) 


Figure 11. Load and frequency deviations (a) random load pattern and (b) frequency deviation in areas with 
optimal ANN controller under random step load change 


6. CONCLUSION 

This study intended to evaluate the performance of optimal ANN controller in a robust LFC problem 
of a four-area interconnected power system. The system is composed with tandem compound steam turbines 
and generation rate constraints under different loading conditions. The Z-N rules are used to find the 
minimum values of ITAE in the design of controllers. The linear quadratic regulator is used to minimize the 
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performance index and hence the mean squared error. The simulation results for the equal change in 
percentage of load with proposed optimal ANN controller gives a significant improvement in terms of time 
response specifications. The settling time, maximum undershoot, and objective function values compared to 
tuned PID and NARMA-L2 controllers. The conventional controllers give higher values of error due to non- 
linearity present in the system and loads. The robustness test of proposed optimal ANN controller is carried 
out with random load pattern. The proposed controller performs satisfactorily under random step load 
changes and thus desirable dynamic control of the system is achieved. 
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